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A MATHEMATICAL  MODEL  OF  PENETRATION  OF  CHUNKY 
PROJECTILES  IN  A GELATIN  TISSUE  SIMULANT 


I.  INTRODUCTION. 

Since  the  late  1930’s  and  early  1940’s  British  and  American  researchers  in  wound 
ballistics  have  been  using  20%  gelatin  gel  as  a tissue  simulant  in  testing  ballistic  projectiles  as  diverse 
as  irregular  grenade  fragments  and  high-velocity  bullets.  Gelatin  is  used  because  it  is  homogeneous, 
presenting  the  same  physical  characteristics  block  after  block;  because  it  is  transparent,  so  that 
events  inside  the  block  can  be  recorded  by  high-speed  movies;  because  its  retarding  properties  are 
similar  to  those  of  skeletal  muscle;  and  because  the  energy  deposit  in  gelatin  correlates  well  to 
measures  of  tissue  damage  and  the  resulting  incapacitation  of  soldiers.  The  disadvantages  of  using 
gelatin  are  that  Firing  tests  are  expensive  and  the  results  are  applicable  only  to  the  weapon  tested. 

A mathematical  model  of  penetration  of  gelatin  tissue  simulant  was  derived  in  a 
previous  report*  and  shown  to  scale  the  penetration  distance  of  a variety  of  spheres  of  different 
sizes  and  densities.  In  this  report,  it  is  shown  that  this  model  may  be  fitted  to  the  data  on 
penetration  versus  time  taken  from  high-speed  movies  of  spheres,  cylinders,  cubes,  and  irregular 
“chunky-shaped”  fragments  penetrating  gelatin.  This  model  makes  possible  an  accurate  prediction 
of  the  energy  deposit,  and  therefore  the  potential  for  incapacitation,  of  projectiles  of  any  size, 
density,  and  striking  velocity,  provided  that  their  shapes  are  similar  to  the  shapes  mentioned  above. 
A brief  review  of  the  derivation  of  the  important  equations  derived  from  that  model  will  be 
presented  in  the  next  section. 

II.  RESULTS. 

The  terms  which  will  be  used  in  developing  the  retardation  models  are  defined  as 

follows' 


Variables: 

F - retarding  force  on  projectile  (dynes) 
t - time  after  impact  (seconds) 
x - distance  penetrated  (centimeters) 
v - velocity  of  projectile  (centimeters  per  second) 

Av  - velocity  loss  at  impact  (centimeters  per  second) 

vQ  - inferred  initial  velocity  in  gelatin  (centimeters  per  second) 


*Sturdivan,  L.  M.  Edgewood  Arsenal  Technical  Report  EB-TR-73022.  A Mathematical  Model  for  Assessing 
Weapons  Effects  from  Gelatin  Penetration  by  Spheres.  September  1973. 
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Constants: 


Gelatin  properties 

b - boundary  layer  thickness  (centimeters) 

M - coefficients  of  velocity  (grams  per  centimeter  second) 
p - density  (grams  per  cubic  centimeter) 

Projectile  properties 

A - mean  presented  area  (square  centimeters) 
m - mass  (grams) 

v§  - striking  velocity  (centimeters  per  second) 
p'  - density  (grams  per  cubic  centimeter) 

Proportionality  (curve  fit) 

a - velocity-loss  coefficient  (grams  per  cubic  centimeter) 
c - velocity-loss  coefficient  (centimeters  per  second) 

Cj  - inertial-force  coefficient  (dimensionless) 

Cv  - viscous-force  coefficient  (dimensionless) 

A The  Retardation  Equation. 


An  application  of  dimensional  analysis  and  elementary  physical  principles  leads  to  the 
proposal  that  the  retarding  force  on  the  projectile  be  considered  the  sum  of  two  components:  an 
inertial  component  which  arises  from  overcoming  the  inertia  of  the  gelatin  which  must  be  moved 
aside  as  the  projectile  penetrates  and  a viscous  component  which  represents  the  friction 
encountered  as  the  projectile  slides  through  the  gelatin.  The  second  component  is  called  viscous 
because  the  gelatin  is  thixotropic;  that  is,  it  liquefies  under  pressure.  Thus,  the  penetrating  projectile 
is  surrounded  by  a boundary  layer  of  viscous  liquid  which  lies  between  it  and  the  solid  gel.  The 
resulting  force  equation  is: 


F = 


dv  Apv  _ n 
_ni__=Cv-— + CjpAv2 


(1) 


where  the  coefficients  Cy  and  Cj  indicate  the  viscous  and  inertial  terms.  The  model  does  not  apply 
at  extremely  high  or  low  velocities  but  it  fits  well  through  a wide  intermediate  range  of  velocities 
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which  is  (he  focus  of  practical  interest.  The  mode!  should  not  be  expected  to  hold  for  penetration 
velocities  approaching  the  speed  of  sound  in  gelatin  (about  1500  m/sec)  since  conipressional  effects, 
which  are  not  modeled,  become  important.  However,  these  transonic  velocities  seldom  occur  in 
practice.  As  the  projectile  nears  the  end  of  its  penetration,  its  veloc.jy  falls  to  a level  so  low  that  the 
pressure  exerted  on  the  gelatin  is  not  enough  to  liquefy  it.  Because  the  projectile  is  then  penetrating 
an  elastic  solid  rather  than  a viscous  liquid,  it  comes  to  a rather  abrupt  stop.  The  errors  in  the  mode! 
at  these  low  velocities  are  safely  ignored  because  of  the  small  amount  of  energy  remaining  in  the 
projectile.  This  model,  a generalization  of  Resal’s  law  (named  after  M.  H.  Resal  who  first  proposed 
this  type  of  force  equation  in  1895),  has  the  following  solution  for  penetration  distance  x as  a 
function  of  time. 


x = 


m 


C,pA 


Cn 


, £l  ^ (i  . 

CV  « \ 


-CyP  At 
bm 


Equation  3 below  gives  velocity  as  a function  of  penetration. 


v = I v 


o+£vm\  /cpA*\  _Cvm 
° Cjpb J \ 1 m ) C,pb 


(2) 


(3) 


where  vc  is  the  inferred  initial  velocity  at  entrance  into  the  gelatin. 

Since  p and  h are  unknown  constants  associated  with  the  gelatin,  they  will  be  grouped 
with  the  Cy  constant  in  the  remainder  of  the  report. 

B.  Impact  Velocity  Loss. 

It  has  long  been  known  that  energy-absorbing  surface  effects,  such  as  baeksplash  and 
the  generation  of  shock  waves  and  surface  waves,  occur  when  a projectile  strikes  the  gelatin 
surface.*  These  effects,  of  course,  are  accompanied  by  a reduction  in  the  projectile  velocity.  It  had 
been  assumed  that,  with  projectiles  as  dense  as  steel,  this  velocity  loss  was  negligible.  However, 
when  equation  2 was  fitted  to  gelatin  time-penetration  data  for  a steel  projectile  with  high  striking 
velocity,  the  resulting  curve  tended  to  over-estimate  the  first  few  points,  suggesting  that  the  initial 
slope  of  the  curve,  which  had  been  presumed  to  be  equal  to  the  striking  velocity,  was  too  high.  In 
other  words,  the  difference,  Av,  between  the  impact  velocity,  vg,  and  the  inferred  initial  velocity, 
vQ,  was  large.  A combination  of  equation  3 and  a similar  equation  for  velocity  as  a function  of  time 
yields  the  following 


[vol  A pA 

— =C,p-X+Cy  — t.  (4) 

v J m bm 

An  iterative  nonlinear  least  squares  scheme  was  used  to  fit  equation  4 to  movie  data  for  the  spheres, 
cubes,  cylinders,  and  fragments  which  are  included  as  the  first  10  projectiles  in  table  1.  This  method 


* McMilien,  Howard  J.  Shock  Wave  Pressures  in  Water  Produced  by  Impact  of  Small  Spheres.  The  Physical 
Review  68,  Numbers  9 and  10  ( 1945). 


9 


Table  1 . Physical  Characteristics  of  Projectiles 


■ 

Projectile 

Materials 

Mass 

Mean 

dimension* 

Mean 

presented 

area 

Density 

gm 

cm 

cm2 

gm/cm3 

0.5-Grain  cylinder 

Steel 

0.0318 

0.175 

0.0361 

7.60 

0.5-Grain  W cylinder 

Tungsten 

0.140 

0.0231 

16.10 

0.85-Grain  sphere 

Steel 

0.055 

0.238 

0.0445 

7.78 

XM36  Fragments 

Steel 

0.065** 

NA 

0.067** 

7.0 

T57  Fragments 

Steel 

0.10** 

NA 

0.087** 

7.0 

2.1 -Grain  cube 

Steel 

0.135 

0.265 

0.1050 

7.31 

16-Grain  sphere 

Steel 

1.041 

0.635 

0.3167 

7.76 

16-Grain  W cube 

Tungsten 

1.020 

0.393 

0.2316 

16.8 

1 6-Grain  cube 

Steel 

1.029 

0.514 

0.3966 

7.57 

225-Grain  cube 

Steel 

14.694 

1.236 

2.2933 

7.77 

7-Grain  sphere 

Steel 

0.439 

0.476 

C.  1781 

7.77 

7-Grain  spheroid 

Tungsten 

0.454 

0.374 

0.1 100 

16.5 

20-Grain  spheroid 

Tungsten 

1.300 

0.545 

0.2334 

15.6 

* Diameter  for  spheres  and  square  cylinders;  edge  for  cubes. 
’’Mean  values. 


fitted  for  v0  as  well  as  the  Resal’s  law  coefficients  Cj  and  Cyn/b.  However,  these  fitted  v0’s  were 
very  poorly  determined  since  v0  is  just  the  slope  of  the  time-penetration  curve  at  zero  penetration 
and  the  curve  is  extrapolated  backward  from  the  data  at  that  point.  A method  was  found  of  pooling 
the  data  from  several  rounds  with  about  the  same  vs  to  determine  a common  vQ/vs  ratio  for  the 
group.*  These  pooled  points  were  then  used  to  derive  a model  and  fit  for  the  required  coefficients. 

Several  suggested  models  of  velocity  loss  due  to  surface  effects  were  found  in  the 
literature  and  examined  for  applicability  to  the  current  problem.  None  were  found  to  be  entirely 
satisfactory,  since  they  either  fitted  the  data  poorly  or  had  improper  boundary  conditions.  Dubin's 
model,**  however,  suggested  a model  of  the  following  type:  Suppose  that  the  momentum  lost  at 


’Details  of  this  methodology  will  be  published  in  a separate  report  entree,  “Consequences  of  Shock  Waves 
Produced  by  Projectile  Impact  on  Tissue." 

’’Dubin,  Henry  C.  Ballistic  Research  Laboratory  Memorandum  Report  2423.  A Cavitational  Model  for  Kinetic 
Energy  Projectiles  Penetrating  Gelatin.  December  1974. 
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impact  is  proportional  to  the  geometric  mean  of  the  impact  and  entrance  energies  and  inversely 
proportional  to  the  density  of  the  projectile;  that  is, 


m(vs  - v0)  = mAv  « ~ 

P V 


2 1 2 
2 mvs  • J mvo 


or 


Av 


v 


o 


This  indicates  that  a plot  of  p'  Av/vQ  versus  v0  would  be  a straight  line.  Instead,  the  plot  shows 
curvature,  suggesting  an  exponential  rise  of  the  form 


(5) 


Solving  for  vQ  in  terms  of  vs,  we  obtain 


v 


o 


The  data  were  fitted  to  equation  6,  yielding  the  values 

a = 0.295  (gm/cm3) 
c = 82,000  (cm/sec). 


(6) 


The  fitted  curve  is  plotted  in  figure  1 together  with  group  mean  values  for  the  supporting  data.  Note 
that  a has  the  dimensions  of  density.  If  we  divide  a by  the  density  of  gelatin  (1.07  gm/cm3  ),  we  get 
a dimensionless  constant  with  a value  0.28.  Note  that  the  model  does  not  distinguish  different  sizes 
or  shapes  of  projectiles  (the  mass  was  in  the  early  equations  but  divided  out  in  the  final 
form  - equation  6).  However,  if  one  examines  the  data  on  individual  rounds  where  the  orientation 
can  be  observed,  as  with  the  large  cubes,  a difference  can  be  seen  between  those  which  struck  nearly 
face-on  and  those  which  struck  more  edge-on.  This  is  because  the  instantaneous  compression,  the 
dominant  feature  of  entry  into  a denser  medium,  is  maximum  when  the  colliding  surfaces  are 
parallel.  Little  use  can  be  made  of  this  fact,  though,  because  all  of  the  projectiles  used  in  this  study 
tend  to  strike  with  random  orientation.  The  best  recourse,  under  these  conditions,  is  to  use  the 


1 1 


mean  curve. 


KEY  TD  PLOT  SYMBOLS 


5/W  -S5D1  A1DG'13A 
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Figure  1 . Velocity  Loss  Data  and  Fitted  Curve  for  Various  Projectiles  at  Impact  on  Gelatin 


C.  Fitting  Che  Ret.  d?tion  Equation. 


The  mcvie  time-- penetration  data  on  the  projectiles  of  table  1 were  again  fitted  to 
equation  4 but  v/ith  vQ  as  a known  parameter  from  equation  6.  The  best-fit  curve  tor  each 
individual  round  was  then  used  to  calculate  velocity  points  corresponding  to  the  time  and 
penetration  points.  These  data  were  then  pooled  with  other  rounds  of  that  projectile  to  fit  for  the 
“pooled”  coefficients  in  equation  4.  This  technique  met  with  a surprising  lack  of  success. 

For  those  projectiles  where  total  penetration  distance  was  known,  the  coefficients  of 
the  median  penetration  round  of  a group  with  relatively  homogeneous  striking  velocity  were  a much 
better  representation  of  the  group  than  the  pooled  coefficients  were.  A select  group  of  these  median 
round  coefficients  was  assembled.  They  seemed  to  fall  into  three  categories:  spheres,  platelet-like 
fragments,  and  everything  else  (cubes,  cylinders,  and  chunky  fragments).  It  is  obvious  from 
equations  2 and  3 that  the  Resal’s  coefficients  are  coupled.  Merely  averaging  the  select  values  of  C’j 
and  Cyp/b  in  the  three  categories  did  not  yield  good  representative  values  of  the  coefficient' 
Several  different  pairs  of  coefficients  from  each  group  were  te.ted  on  the  rest  of  the  group.  Those 
that  did  best  in  the  entire  select  group  of  median  penetration  rounds  were  clustered  about  the 
rounded  off  values  given  in  table  2.  These  values  were  then  tested  against  time-penetration  uiia 
including  rounds  with  penetration  higher  and  lower  than  the  median.  As  expected,  the  predicted 
curves  lay  below  or  above  the  data  in  those  cases.  However,  approximately  equal  numbers  fell  on 
either  side  for  each  projectile.  As  mentioned  earlier,  that  portion  of  the  data  where  the  projectile 
did  not  liquefy  the  gelatin  and  abruptly  stopped  was  not  considered  in  judging  the  fit  of  the  curves. 
This  generally  occurs  at  a velocity  between  50  and  100  m/sec. 


Table  2.  Coefficients  for  the  Resal’s  Law  Equation 


Projectile  type 

Cl 

Cyp/b 

Spheres 

0.10 

3000 

Cubes,  cylinders,  fragments 

0.175 

3000 

Platelets  (XM36  only) 

0.15 

5000 

Figures  2 through  4 show  some  of  the  exact  fits  to  Resal’s  law.  In  these,  the  values  of 
C|  and  Cyp/b  are  unique  for  each  round.  Figures  5 through  14  show  data  from  median  penetration 
rounds  plotted  against  curves  using  the  general  coefficients  from  table  2.  Figures  15  and  16  show 
predicted  general  curves  versus  data  from  a few  high-  and  low-penetration  rounds.  Figures  1 7 
through  19  show  predicted  curves  versus  data  for  a group  of  spheres  and  spheroids  not  used  to 
derive  sphere  coefficients  but  used  to  test  them  before  inclusion  in  table  2. 
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Figure  3.  Resal's  Law  Curves  Individually  Fitted  to  Penetration-  Time  Data  for  Steel  Cubes 
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Figure  6.  Generalized  Curve  Versus  Penetration-Time  Data  for  the  0.5-Grain  Tungsten  Cylinder 
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Figure  7.  Generalized  Curve  Versus  Penetration- Time  Data  for  the  0.85-Grain  Steel  Sphere 
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Figure  8.  Generalized  Curve  Versus  Penetration-Time  Data  for  the  2.1 -Grain  Steel  Cube 


Figure  1 1.  Generalized  Curve  Versus  Penetration-Time  Data  for  the  16-Grain  Tungsten  Cube 
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Figure  1 2.  Generalized  Curve  Versus  Penetration-Time  Data  for  the  16-Grain  (1/4-Inch)  Steel  Sphere 


Figure  13.  Generalized  Curve  Versus  Penetration-Time  Data  fot  the  XM36  Platelet  (Fragment) 
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Figure  14.  Generalized  Curve  Versus  Penetration-Time  Data  for  the  T57  Preformed  Fragment 
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Figure  15.  Generalized  Curve  Versus  Penetration-Time  Data  for  Several  High-Penetration  Rounds 
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Figure  17.  Generalized  Curve  Versus  Penetration-Time  Data  for  the  7-Grain  (3/16-lnch)  Steel  Sphere 


Figure  19.  Generalized  Curve  Versus  Penetration-Time  Data  for  the  ^0-Grain  Tungsten  Spheroid 


Figures  20  through  24  show  some  time- penetration  data  from  irregular  cast  iron  shell 
fragments.  Before  they  were  fired,  these  fragments  were  visually  sorted  into  several  categories  of 
shapes.  Figures  20  through  22  represent  the  “chunky”  or  compact  fragments,  whereas  figure  23 
represents  the  “very  irregular”  category  composed  of  fragments  with  a very  irregular  surface  having 
lumpy  projections  randomly  extruding  outward  from  it  and  figure  24  represents  the  long 
“splinterlike”  fragments  often  seen  recovered  from  exploded  cast  iron  shell.  These  data  are 
contrasted  with  the  curve  predicted  for  fragments  of  like  mass,  presented  area,  and  velocity  with  the 
generalized  cube/fragment  coefficients  from  table  2.  Physical  characteristics  of  these  projectiles  are 
listed  in  table  3.  Note  that  the  mass  and  mean  presented  area  of  each  fragment  are  unique.*  Because 
it  is  assumed  in  the  model  that  each  projectile  has  constant  mass  and  mean  presented  area,  data  on 
fragments  which  broke  into  two  or  more  pieces  upon  impacting  the  gelatin  were  not  used  in  the 
analysis  nor  in  the  figures. 


Table  3.  Characteristics  of  Irregular  Cast  Iron  Fragments 


Fragment 

No. 

Category* 

Mass 

Mean  presented 

area 

gm 

cm2 

5 

Very  irregular 

1.23 

0.471 

1 i 

Very  irregular 

0.95 

0.419 

16 

Very  irregular 

0.787 

0.374 

19 

Very  irregular 

0.439 

0.265 

27 

Chunky 

4.55 

1.155 

28 

Chunky 

4.83 

1.077 

29 

Chunky 

3.48 

0.936 

30 

Chunky 

2.94 

0.761 

31 

Chunky 

3.06 

0.794 

32 

Chunky 

2.89 

0.794 

37 

Chunky 

2.33 

0.652 

40 

Chunky 

2.62 

0.652 

45 

Chunky 

1.31 

0.452 

48 

Chunky 

0.683 

0.290 

49 

Chunky 

0.793 

0.346 

51 

Chunky 

0.652 

0.265 

55 

Chunky 

0.728 

0.290 

58 

Chunky 

0.283 

0.136 

64 

Long 

4.827 

1.284 

68 

Long 

1.349 

0.561 

71 

Long 

0.485 

0.342 

95 

Chunky 

0.170 

0. 1 85 

*See  text. 


*Mean  presented  areas  of  the  fragments  were  measured  on  the  automatic  shell-fragment  area-measuring  device 
(ASFAM-D)  at  the  Materiel  Test  Directorate,  Aberdeen  Proving  Ground. 
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VERY  IRREEULRR  FRREMENT5 


Figure  23.  Generalized  Curve  Versus  Penetration-Time  Data  for  Several  “Very  Irregular”  Fragments 


Figure  24.  Generalized  Curve  Versus  Penetration- Time  Data  for  Several  “Splintertike”  Fragments 


III.  CONCLUSIONS. 


The  velocity  loss  and  Resal’s  law  coefficients  were  derived  on  the  basis  of  abundant 
data  on  steel  and  tungsten  projectiles  and  were  shown  to  be  reasonably  good  predictors  of 
penetration  (except  for  the  different  phase  portion  near  the  stopping  point)  for  a wide  range  of 
shapes  and  over  two  orders  of  magnitude  in  mass.  Particularly  gratifying  is  the  ability  of  the  model 
to  predict  the  very  slight  difference  in  penetration  of  the  small  steel  cylinders  at  i500-  versus 
20C0-m/sec  striking  velocity  (see  figure  5).  This  phenomenon  was  previously  considered  an  anomaly 
introduced  by  the  much  greater  deformation  caused  at  the  higher  velocity  impacts.  Although  that 
deformation  probably  does  account  for  the  more  rapid  halt  of  the  2000-m/sec  rounds  (shown  by  an 
earlier  deviation  from  the  predicted  curve  and  shorter  overall  penetration),  it  is  seen  that  moderate 
deformation,  without  breakup,  does  not  cause  such  a deviation  from  the  model  that  its  usefulness  is 
lost. 


Some  caution  should  be  taken  in  using  these  models  to  extrapolate  far  beyond  the 
range  of  physical  characteristics  of  the  projectiles  from  which  the  models  and  constant  parameters 
were  derived  and  especially  *:o  very  high  velocity  impacts.  Compression,  deformation,  and  breakup 
increase  rapidly  with  increasing  velocity  and  can  completely  invalidate  the  model.  Within  these 
constraints,  the  model  may  be  used  to  calculate  functions  of  penetration  distance  or  time,  such  as 
velocity,  energy  deposit  and  acceleration,  to  be  used  in  weapons  assessments  or  war  game  models. 

Although  no  data  were  available  at  the  time  of  this  writing  on  time-penetration  into 
gelatin  by  the  less  dense  projectiles,  the  velocity  loss  inferred  by  Dubin*  from  gelatin  cavity 
measures  on  penetrating  nylon  and  aluminum  spheres  agrees  very  well  with  predictions  by  the 
present  model.  As  Dubin  assumed  the  same  drag  coefficient  for  these  data  as  he  did  for  those  of 
steel,  we  interpret  those  results  as  generally  favorable  for  the  present  model. 

It  is  not  anticipated  that  low-density  materials  will  be  extensively  used  in  future 
weapons.  However,  the  increasing  use  of  light  materials  in  vehicles,  particularly  the  use  of 
lightweight,  high-strength  alloys  in  armoring  those  vehicles,  will  make  spall  injuries  by  low-density 
fragments  much  more  common  on  future  battlefields.  This  suggests  that  limited  firings  of 
low-density  fragments  would  be  worthwhile  to  test  the  accuracy  of  these  models  in  predicting  their 
penetration  and  energy  deposit. 


*Dubin,  Henry  C.  Ballistic  Research  Laboratory  Memorandum  Report  2423.  A Cavitational  Model  for  Kinetic 
Energy  Projectiles  Penetrating  Gelatin.  December  1974. 
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